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A hallmark feature of topological physics is the presence of one-way propagating chiral modes at 
the system boundary [1, 2]. The chirality of edge modes is a consequence of topological character 
of the hulk. For example, in a non-interacting quantum Hall (QH) model, edge modes manifest as 
mid-gap states between two topologically distinct bulk bands. The bulk-boundary correspondence 
dictates that the number of chiral edge modes, a topological invariant called the winding number, is 
completely determined by the bulk topological invariant, the Chern number [3]. Here, for the first 
time, we measure the winding number in a two-dimensional (2D) photonic system. By inserting a unit 
flux quantum to the edge, we show that the edge spectrum resonances shift by the winding number. 
This experiment provides a new approach for unambiguous measurement of topological invariants, 
independent of the microscopic details, and could possibly be extended to probe strongly correlated 
topological orders. 

Recently, there has been a surge of interest in investigating topological states with synthetic gauge fields. 
Synthetic gauge fields have been realized in various atomic [4—7] and photonic systems [8, 9]. In particu¬ 
lar, topological photonic edge states have been imaged in two recent concurrent experiments [10, 11] and 
robustness of their transport has been quantitatively confirmed both in the microwave [12] and telecom 
domains [13]. Several other interesting proposals have investigated topological states in ID [14-16], 2D 
[17-21] and also 3D [22] synthetic structures. Topological states are characterized by topologically invari¬ 
ant integers [2]. More specifically, the bulk-boundary correspondence dictates that the number of chiral 


edge modes, the winding number, is completely determined by the bulk topological invariant, the Chern 
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number [3]. In fermionic systems, conductance measurements reveal these quantum numbers. However, 
a direct measurement of these integers is non-trivial in bosonic systems, mainly because the concept of 
conductance is not well defined [23, 24]. While these integers have been measured in ID bosonic systems 
[15, 25, 26], the 2D bosonic case has been limited to atomics lattices [7]. 

Here, we experimentally demonstrate that selective manipulation of the edge can be exploited to measure 
topological invariants, i.e., the winding number of the edge states. We implement an integer quantum Hall 
system using a fixed, uniform synthetic gauge field in the bulk and couple an additional tunable gauge field 
only to the edge. The edge state energy spectra flows as a function of this tunable flux. With insertion of 
a unit quantum of flux, the edge state resonances move by ±1, which is the winding number of edge states 
in our system. This spectral flow can be directly observed in an experiment as the flow of transmission 
resonances, and thus provides a direct measurement of the winding number of edge states. For this demon¬ 
stration, we employ the unique ability of our photonic system to selectively manipulate edge states — a 
feature that is challenging to achieve in current electronic and atomic systems. 

To model the spectral flow of a quantum Hall edge, with winding number k = 1, we consider a linear 
edge dispersion Ep = vp where Ep is the energy, v is the group velocity, and p is the momentum along the 
edge. When a gauge flux (6) is coupled to the edge, the momentum is replaced by the covariant momentum, 
i.e.. 


Ep 



( 1 ) 


where L is the length of the edge and q is the charge of the edge excitations. For non-interacting photons, the 
charge <7 = 1. Note that the corresponding vector potential is simply 0/L. For a finite system, quantization 
of momentum on the edge results in 
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where n is an integer. Thus, insertion of 0 = 27r flux shifts E^ —> E^-i resulting in an anomalous spectral 
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FIG. 1: a, Spectral flow of edge states with coupled gauge flux 0 for g = 1 and k — 1. When 9 = 2tt, edge states 
move by fc = 1 positions. b,c, Square annulus of coupled ring resonators which implements the integer quantum 
Hall effect, with outer edge of 10 sites and inner edge of 4 sites. The uniform synthetic magnetic field for photons 
is introduced by asymmetrically positioning the link rings. The system exhibits two pairs of clockwise and counter 
clockwise propagating edge states, one each on the outer and the inner edges. The tunable gauge field coupled only 
to the edge states is introduced by fabricating heaters on link resonators situated on the lattice edges. The total flux 9 
introduced by heaters is the sum of individual phases incurred in link resonator arms, d, SEM image showing heaters 
fabricated on top of link resonators (green), e. Schematic of the waveguide cross section showing the ring resonators, 
the metal heaters and the metal routing layer. 


flow, as shown in Fig. la. This is in contrast to the case of a topologically trivial system, where there is 
no net shift and having a zero and 27r fluxes are equivalent to each other. However, in QH systems, these 
linear-dispersion mid-gap modes continuously interpolate between two topologically distinct bulk bands 
and the difference between the Chern numbers of these bulk bands dictates the shift, and hence the winding 
number of the edge states. For winding numbers larger than one (k > 1), there exist k edge modes in 
the bandgap and hence k copies of this relation, i.e., for a unit quantum of flux insertion, the edge state 
resonances shift by k units. Such spectral flow is similar to Laughlin’s charge pump [27, 28] but with an 
important distinction: here, the gauge flux is coupled only to the edge of the system (see SI). Moreover, 
the observable in our interferometric measurement is expected to be an integer, in contrast to the scheme 
presented in [23] which measures Hall drift, a continuous variable. More generally, our scheme provides a 
powerful universal probe (independent of microscopic details) of topological order that can be generalized 
to the case of strongly correlated topological systems. 

To experimentally observe and measure this spectral flow, we implement the integer quantum Hall model 
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FIG. 2: a, Measured and b, simulated transmission (T, in linear scale) as a function of phase 9 and ut (frequency), for 
spin-up excitation. For a 27r increase in gauge flux, the edge state resonances move by one position, giving k = -|-1.0(1) 
for CW edge band and k = —1.0(2) for the CCW band. Insets: Zoom-in of the edge state bands, c. Measured and d, 
simulated transmission spectra (log scale) for 6/2tt — 0,1. For this 27r increase in flux, the measured spectra match 
approximately in the edge state regions. The green and red shaded regions indicate the CW and CCW edge state bands, 
respectively. The transmission is normalized such that the maximum is unity. Here ujq is the resonance frequency of 
the resonators and J is the coupling rate between site resonators. 


in a photonic system: a two dimensional square annulus of ring resonators with a uniform synthetic mag¬ 
netic field in the bulk and a tunable gauge field coupled only to the lattice edge (Fig. Ib-e). The uniform 
magnetic field with flux (f)Q = '^ radians per plaquette is synthesized using asymmetric placement of site 
and link resonators, as previously described in Ref. [10]. To couple a tunable gauge field to the edges, we 
fabricate metal heaters above the link ring waveguides on the lattice edge (Fig. Ic-e). These heaters use 
thermo-optic effect to modify the accumulated phase of light propagating through the waveguides and hence 
result in a gauge flux. 

Our system supports two pseudo-spin components, up and down, which circulate in opposite directions 
in ring resonators [10]. Each pseudo-spin component experiences a uniform magnetic held, and we can 
selectively excite each pseudo-spin by appropriately choosing the input port (Fig. Ic). Furthermore, with 
this particular choice of uniform magnetic held, the energy spectrum exhibits two bandgap regions [3]. For 
spin-up excitation, the low energy bandgap corresponds to clockwise (CW) circulating edge states confined 
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on the outer edge, and counter-clockwise (CCW) circulating edge states on the inner edge, as shown with 
green arrows in Fig. lb with their corresponding winding numbers. Similarly, the high energy bandgap 
corresponds to edge states both on the outer and the inner edges, however, they circulate opposite to those 
in the low energy bandgap, as shown in Fig. lb by red arrows. A pseudo-spin flip swaps the position of 
CW and CCW edge bands in the transmission spectrum, and also flips the sign of gauge flux 9 [29]. 

Fig. 2a,b show the measured and simulated transmission spectra as a function of the coupled flux 6, for 
spin-up excitation. Edge states of the outer edge (bright regions) and the bulk states (dark region in the 
middle) are easily identifiable. As the coupled flux 9 increases, the energy of the CW edge states decreases, 
whereas the energy of CCW edge states increases. For a 27r increase in flux, the edge state resonances move 
by one resonance to replace the position once held by its neighbor. This flow indicates that the measured 
winding number is fc = -|-1.0(1) for the CW circulating edge states, and k = —1.0(2) for the CCW 
circulating edge states. We note that only the outer edge contributes to the transmission spectra, because the 
edge states are tightly confined to the lattice boundary. Moreover, such sharpness of edge states allows us to 
locally create the extra gauge flux 9 by placing the heaters only on the outmost rings at the boundary. Also 
note that the flow of resonances observed here is very distinct in appearance and physical origin from what 
is observed in a Fabry-Perot cavity, where increasing the length of the cavity always reduces the resonant 
frequencies. However, in our observation CW and CCW edge state resonances shift in opposite directions, 
and resonance peaks shift exactly by one. 

Fig. 2c,d show the overlap of the observed and simulated transmission spectra at ^ = 0 and 1, show¬ 
ing excellent qualitative agreement with the simulation results. We attribute the small discrepancy in the 
overlap to lattice disorder introduced by the nonlocal response of the heaters (see SI). The discrepancy is 
more pronounced for the CCW state because of finite size effects in travelling a shorter path from input to 
the output port. Also note that compared to measured spectra, simulated spectra show high transmission 
through the two bulk bands. This is because the bulk bands are not topologically protected and therefore, 
disorder in the experimental system suppresses transmission through these bands. 

In order to verify the local character of the gauge coupling, we selectively couple a nonzero flux 9 to the 
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FIG. 3: a, Measured and b, simulated transmission spectrum (linear scale, normalized to unity) as a function of flux 6 
when only the inner edge is heated. The energies of the edge states associated with the outer edge do not change, c, 
Measured and d simulated transmission spectrum, for spin-down excitation. The CW and CCW circulating states have 
now exchanged positions and move in opposite directions (compared to Fig. 2a,b), giving k = 1.2(2) for CW edge 
band and k = —1.2(2) for the CCW band. The transmission spectra shown in a and c are from two different devices. 


inner edge by heating the link rings on the inner edge. We observe that heating only the inner edge does 
not shift the edge states (Fig. 3a,b), which indicates the strong confinement of edge states at the lattice 
boundary. Furthermore, we excited the system with a flipped pseudo-spin (compared to Fig. 2). Fig. 3c,d 
shows the measured and simulated transmission spectra as a function of coupled flux 9. Since the spin flip 
reverses both the sign of the coupled flux 6 and the position of edge bands in the transmission spectrum, the 
resulting spectrum is similar to Fig.2. 

As a control experiment, to investigate the difference of gauge coupling to chiral and non-chiral states, 
we fabricate a ring of resonators (Fig. 4a) which corresponds to a ID tight-binding model. This system 
supports topologically trivial Bloch modes, circulating in CW and CCW directions, which are degenerate in 
the absence of a magnetic flux. Fig. 4b,c show the experimentally observed and the simulated transmission 
spectra as a function of coupled flux 6. The flux lifts the degeneracy of the CW and CCW states and with 
an increase in the flux, these states move in opposite directions. Therefore, for 9 — 2tt, there is no net 
shift of the resonances, with measured k = 0.03(8). The absence of spectral flow proves that this system is 
topologically trivial. As shown in Fig. 4d, we observe the transmission spectra at 9/2Tr = 0,1 exhibit an 
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FIG. 4: a A ring geometry to investigate gauge coupling to topologically trivial states. The ring supports CW and 
CCW circulating modes. External flux can be coupled to these states by using heaters on the link resonators, b, 
Measured and textbfc, simulated transmission as a function of external flux 9 and frequency ui. The CW and CCW 
modes move in opposite directions but for a unit increase in flux, they do not change position, d, Measured and e, 
simulated transmission spectra at 0 = 0, 27r overlap well. 


excellent overlap. However, compared to the simulation, we realize that not all the states can be resolved in 
the experimental data. Such discrepancy is due to the absence of topological protection of non-chiral states 
against disorder which leads to coupling between CW and CCW modes. 

Here, we used a simple model to describe spectral flow of a chiral edge coupled to a gauge field, for 
non-interacting bosons. However, the emergence of anomalous spectral flow is more general [1,2], and can 
be modeled using the elegant framework of conformal field theory of chiral bosons, where this behavior 
is known as chiral anomaly [1, 30, 31]. This framework allows the investigation of systems with strong 
photon-photon interaction, such as circuit-QED [32] and polariton [33] systems, and also ultracold atoms 
with topological features [4—6, 25]. In particular, the presence of strong interactions could lead to the 
emergence of various fractional quantum states in the steady-state regime [34]. Spectral flow could be an 
interesting way to explore the fascinating many-body features of topological states in systems with synthetic 


gauge fields. 






















Methods 


Heaters: The heaters used in these devices are 110 nm thick metal (Ti) pads, with typical resistance of 
117 Ohms, such that a current flowing through them generates heat and modifies the refractive index, and 
the accumulated phase of light propagating through the waveguide. A typical heater requires « 76.5 mW 
to introduce a 2tt phase. Instead of employing a single heater which may introduce a significant disorder in 
the lattice, we distribute the heaters along the edges while maintaining an accumulated flux of 2tt (see SI). 

Measurement of winding number: To estimate the value of winding number k, we measure the fre¬ 
quency shift of an edge state resonance for 2tt flux insertion; and its separation from the neighboring 
resonance at 6* = 0. The integer k, in a given band, is then the mean of the ratio of shift of resonances to 
their separation. The standard deviation of this ratio is the error on k. 
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I. COUPLING EXTERNAL GAUGE FLUX 

The implementation of a uniform synthetic magnetic field for photons and coupling an external gauge 
field to the edge states hinges on the use of link resonators to couple the site resonators. By vertically 
displacing the link resonator we can introduce a direction dependent, fixed, hopping phase for the photons 
hopping between site resonators [1, 2]. The fabrication of heaters above the link resonator waveguides 
enables us to tune this hopping phase. Here, we analyze a system of two site resonators coupled by a link 
resonator, with a heater on the link resonator waveguide. The heater introduces a phase 29. We show this 
system is equivalent to two site rings coupled by an effective coupling rate J and a hopping phase 9t, where 
9t includes the fixed phase incurred due vertical shift of the link resonator and the tunable external gauge 
flux introduced by the heater. 

We first analyze the system of site rings coupled by a link ring, using the rigorous transfer matrix method. 
The site resonators have a length L and the link ring L + rj. Eurthermore, the link ring is vertically displaced 
by length ^ and has a heater on the upper arm which introduces a phase 26. We label the fields as shown in 



FIG. SI: (a) A system of two site resonators, coupled by a link resonator with a heater, and also coupled to input and 
output waveguides. The figure labels the fields for the transfer matrix analysis, (b) Effective description of the system 
in coupled mode theory. The effect of the link resonator is contained in the coupling rate J and hopping phase 9t which 
accounts for the vertical shift of the link ring and the phase introduced by the heater. 
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Fig. SI. Using transfer matrix formalism, we can relate the different field amplitudes as [2] 
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Here, ti and ki are the field transmission and cross coupling coefficients for coupling of link rings to the 
input and output ports, and t, k are those for the coupling between the rings. /3 and a are the propagation 
constant and the loss coefficient in the ring resonator waveguide, respectively. These equations then give 
the field at the drop port as 


Sd = 




^g*2egi/3(L+,)) (goL^ _ gi/SL^.) + (1 - ft,) 


(S8) 

Here we have assumed that the extra length rj of the link ring and its vertical shift ^ are negligible compared 
to L, so that the loss incurred in these extra lengths is insignificant. In our experimental system, p = 320 
nm, ^ = 80 nm and L « 70/im, which justifies the above assumption. 

In the weak coupling limit, i.e., when Ki,K 1 and ti,t ~ 1, the expression for the drop field can be 
simplified to 
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This system of two site rings coupled by a link ring can be equivalently described by a system involving 
only the site rings, coupled by coupling rate J and a hopping phase zL6t (Fig. Sib). Specifically, we use 
single mode approximation for the site rings and employ the coupled mode theory to capture the effect of 
the link ring in J and 9t- In the coupled mode theory, the rate equations for the time evolution of the site 
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ring energy amplitudes, ai{t) and a 2 {t), are 
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where Hex is the coupling of the site rings to the input and output waveguides and Km is the resonator loss 
rate. A steady-state solution of the above equations for a plane wave excitation of frequency w gives 
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Then, the field output at the drop port = y/2KexCi2 is 
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To compare this expression for the drop field to that derived using the transfer matrix approach, we use 
the relations [2] 
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The drop field in the transfer matrix formulation is then 
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where 9t = 2/3^ + 6. We see that the introduction of a phase 26 in one of the arms of the link resonator 
results in an additional hopping phase of 9. Further, for [dr] + 29 = (tt, Stt, Btt, ...), i.e., when the link ring 
is anti-resonant to the site rings, the two expressions are identical. Because the link ring is anti-resonant, it 
does not store any energy and simply acts as a waveguide, without affecting the coupling between the site 
rings. 

When /3y -f 29 is not an odd-integer multiple of tt, the two expressions for the drop field, derived using 
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FIG. S2: (a) Shifting of eigenvalues as a function of external magnetic flux 8. The annulus now supports two sets 
of edge states, at the outer (b,e) and the inner edges (c,d). In a given bandgap, the outer (f,i) and the inner (g,h) 
edge states circulate around the lattice in opposite directions because of their opposite group velocity. The arrows 
indicate direction of group velocity (energy flow). Furthermore, in the energy spectrum (a), the clockwise and the 
counterclockwise circulating states move in opposite directions because of their winding number which is ± 1. 


the coupled mode theory and the transfer matrix method, are identical if we use an effective coupling rate 
= J! sin and shift the resonance frequency of the site rings as Wg® = wq + J cot . in 

our experiment, we design the link rings such that firj = tt. Further, we distribute the external gauge flux 9 
over 12 heaters so that the introduction of external gauge flux 9 results in negligible correction to coupling 
rate J and a small shift in the resonance frequency of the site rings. For example, for 9 — 2tt divided over 
12 link rings, the correction to the coupling rate J is only. 3.5%. However, the shift in resonance frequency 
Awg is 0.27J, which adds to disorder in the lattice. 

Our system supports two degenerate spin states - circulating CW and CCW in the ring resonators. The 
above discussion is for CCW state. Following similar procedure for CW state (spin flipped) results in a 
negative sign for the gauge flux theta. Also note that coupling a gauge flux by locally heating an arm of 
the link resonator is not the same as increasing the length of the resonator. In particular, the gauge flux 9 
introduced here has opposite sign for opposite direction of hopping (right or left, see Fig.SI b). 
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II. EIGENVALUES OE SQUARE ANNULUS AS A EUNCTION OE EXTERNAL FLUX 

The 2D annulus geometry of site resonators (shown in Fig. 1) can be described using coupled-mode 
theory with a Hamiltonian [1,2] given as 


iTo = - J V 


n p '^y4^0 ,y;edge _|_ At A ^ 

/ ^ ^ ' ^x,y^x-\-l,y^ 


oT'^4>tSx ,y\edge 


x,y 




(520) 

(521) 


where ax,y and dj. y are the photon annihilation and creation operators, respectively, at a site resonator with 
position index x, y, and J is the hopping rate for photons. The hopping phase (f)Q along the x—axis gives 
rise to a uniform magnetic field in the bulk with flux (pQ per plaquette [1,2]. In addition to this uniform 
bulk field, we couple a tunable gauge flux 0 only to the outer (or the inner) edge of the system (Fig. 1). 
The tunable hopping phase (pt = ■^, where iV = 12 is the number of heaters used on the edge to achieve 
an accumulated flux 9. The sign of the tunable phase term is —ve for top and left edges of the annulus; and 
+ve for bottom and right edges of the annulus. 

Fig. S2 shows the eigenvalues of this annulus as a function of the external flux 9. This system supports 
two sets of CW and CCW propagating edge states, one each on the outer and the inner edges. Furthermore, 
edge states on the outer edge (e.g. CW) are paired with edge states of opposite chirality (CW) on the 
inner edge, i.e., they occupy the same bandgap (left). To illustrate the spectral flow of these edge states 
as a function of external gauge flux, we couple the gauge flux to both the outer and the inner edges of 
the system. We observe that with an increase in the external flux, all the edge states shift outwards or 
inwards. The CW-propagating outer edge states in the first bandgap move toward lower energies, whereas 
the CCW-propagating outer edge states in the third bandgap move toward higher energies. The inner edge 
states move exactly opposite to the outer edge states. More importantly, for a 27r increase in magnetic flux, 
each set of states shift exactly by one, their winding number. The direction of shift indicates the sign of 
the winding number, as is shown in Fig. S2. After a 27r increase in the flux, the spectrum returns back 
to its original shape with 9 = 0. Note that in the measured transmission data shown in the main text, we 
observe transmission resonances corresponding only to the outer edge states. This is because in our system, 
the input and output waveguide couplers are coupled only to the outer edge and their coupling to inner 
edge is exponentially suppressed (Fig. 1). Furthermore, it is instructive to compare the flow of resonances 
observed here to that in a Fabry-Perot cavity. As we discussed, introduction of an additional gauge flux is 
not the same as increasing the length of the resonator. In particular, the sign of the gauge flux is opposite for 
CW and CCW propagating states, and therefore, these resonances move in opposite directions. In contrast, 
increasing the length of the Fabry-Perot cavity shifts all resonances towards lower frequencies. 

Also note that in the transmission spectra presented in the main text (Fig. 2, 3), only the outer edge 
contributes to the transmission spectra. This is because the coupling between the probe waveguides to the 
inner edge states is exponentially suppressed as where d is the number of rings between the outer 

and the inner edge (in this case 4) and (g = 1/v^^ ~ 0-8 is the magnetic length in units of the lattice 
spacing. 


III. HEATER CALIBRATION 


We used an add-drop filter (ADF - a ring coupled to two waveguides) with a heater, to calibrate the phase 
shift acquired for a unit heater power. The ADF ring has exactly the same dimensions as the rings in the 2D 
array. The ring waveguides are designed to be 510 nm wide to ensure single mode (TE) operation at telecom 
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Heater Power (mW) 


FIG. S3: Measured phase shift in an ADF, as a function of the heater power. The acquired phase shift (in units of Itt) 
increases linearly with the heater power (in mW), with a slope of 0.026(1). 


wavelengths and the typical ring perimeter is « 70/rm. The coupling gap between the ring waveguides is 
140 nm which gives a coupling rate J = 40.6 GHz. The heater introduces a phase shift 9 in the ADF ring 
waveguide and shifts its resonance frequency by Awq, such that 


9 = 


Aotp 

FSR 


27r, 


(S22) 


where FSR = 27r^ is the free spectral range of the ring, Vg is the group velocity and L is the length of 
the ring. Measuring the resonance frequency shift as a function of heater power gives us an estimate of the 
phase difference introduced by the heater. We used an ADF ring which had exactly the same dimensions 
as the rings in the square annulus and the chain of resonators. Fig. S3 shows the observed phase shift as 
a function of heater power. The heater resistance was «:! 120 Ohms and the voltage applied to the heater 
ranged from 0-3 Volts. We see that the acquired phase increases linearly with the heater power. 

For the square annulus and the chain geometry, we used this calibration curve to estimate 9 for a given 
voltage applied to the heaters. For the annulus, all the 12 heaters in the outer (and the inner) loop (Fig. Ic)) 
were serially connected. The measured resistance was 1.4 KOhms and the applied heater voltage across the 
loop ranged from « 8-14 Volts, which introduced the phase 9 = n — Stt. Similarly, in the chain of rings 
(Fig. 4a), all the 20 heaters were connected serially, giving a total resistance of 2.31 KOhms. The total 
voltage applied across the heaters ranged from « 0 — 13 Volts, for 0 = 0 — 27r. Furthermore, as we saw in 
the previous section, the transfer matrix analysis of a system of two site rings, coupled by a link ring, shows 
that if the heater introduces an extra phase 29 in one of the arms, the effective hopping phase (using CMT) 
between the link rings is 9. We have included this factor of 2 in our calculations. 


IV. DISORDER ESTIMATION 

In our design, the heaters are placed 600 nm above the link resonator waveguides. Because of this, the 
heat is not conhned only to the ring beneath the heater but instead also heats the neighboring site resonators. 
Using a three ring device where two site rings are coupled by a link ring, we estimate that for every unit shift 
in resonance frequency of the link ring, the site rings move by 0.3 units. The small curvature observed in the 
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shifting of edge states is also a result of this disorder. We have accounted for this disorder in our simulation. 
Moreover, because of this disorder we find that for the square annulus, the transmission spectrum retains its 
shape in the range 0 « tt — Stt. In our measurements reported in Fig. 2, the external gauge flux has been 
normalized to 6* = 0 — 27r. 

Furthermore, because of this nonlocal behavior of the heaters, the device spectrum shifts towards lower 
frequencies as we increase the voltage applied to the heaters. We have adjusted for this overall shift of the 
transmission spectra. For each spectrum, we choose two edge state peaks, one each on CW and CCW bands 
and align all the spectra to the center of these two peaks. The resulting plot thus shows only the differential 
shift between the two edge state bands. For the ring geometry, the measured spectrum at 0 = 27r in Fig. 
4d has been scaled along the frequency axis (by a factor of 0.98), to offset the dispersion effect (due to the 
overall shift) which results in slight broadening of the spectrum. 


V. COMPARISON WITH LAUGHLIN’S CHARGE PUMP 

The spectral flow observed here is very similar to Laughlin’s charge pump [3, 4]. The Laughlin’s pump 
was experimentally observed for ID dimensional photonic systems [5, 6], and proposed as a way to measure 
winding number in 2D photonic systems [7]. However, our scheme requires manipulation and probe only 
at the boundary of the system, unlike Laughlin’s case where the gauge is coupled to both the edge and the 
bulk. Therefore, our scheme avoids technical complications that may arise due to the bulk manipulation 
which is required in Ref. [7]. 
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